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Semiclassical theory of proximity effect predicts a gap Eg ~ hD/L 2 in the excitation spectrum of a 
long diffusive superconductor/normal - metal/superconductor (SNS) junction. Mesoscopic fluctuations lead 
to anomalously localized states in the normal part of the junction. As a result, a non-zero, yet exponentially 
small, density of states (DOS) appears at energies below E g . In the framework of the supermatrix nonlinear 
cr-model these prelocalized states are due to instanton configurations with broken supersymmetry. The exact 
result for the DOS near the semiclassical threshold is found provided the dimensionless conductance of the 
normal part Gn is large. The case of poorly transparent interfaces between the normal and superconductive 
regions is also considered. In this limit the total number of the subgap states may be large. 
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PACS: 73.21.-b, 74.50.+r, 74.80.Fp 

1. Introduction. It has recently been shown 
within several different although related contexts that 
the excitation energy spectrum of Superconductive - 
Normal (SN) chaotic hybrid structures pi J3] and su- 
perconductors with magnetic impurities pi 0] does not 
possess a hard gap as predicted by a number of pa- 
pers J|, [| 0, H using the semiclassical theory of su- 
perconductivity §0 0. With mesoscopic fluctua- 
tions taken into account, the phenomenon of soft gap 
appears: the density of states is nonzero at all energies, 
but it decreases exponentially fast below the semiclas- 
sical threshold E g ~ S/ r o with t c being the charac- 
teristic dwell time in the N region. In particular, for 
diffusive systems perfectly connected to a superconduc- 
tor, Eg has the order of the Thouless energy i?Th in the 
N region § §. 

The first result in this direction was obtained in 
Ref. Q], where the subgap DOS in a quantum dot 
was studied by employing the universality hypothesis 
and the predictions jL2| of the random matrix the- 
ory (RMT) Later on, the tail states in a super- 
conductor with magnetic impurities were analysed in 
Refs. Q on the basis of the supersymmetric nonlinear 
cr-model method fl4|| extended to include superconduc- 
tive paring |Te{ ]. 

Fully microscopic approach to the problem of the 
subgap states in diffusive NS systems was developed in 
Ref. in the framework of the supersymmetric cr-model 
similar to that employed in Refs. (|[ [|. Physically, the 
low-lying excitations in SN structures are due to anoma- 
lously localized eigenstates |ll| in the N region. From 
the mathematical side, nonzero DOS at E < E g comes 
about when nontrivial field configurations - instantons 



- are taken into account in the cr-model functional inte- 
gral. As shown in [0 , at E w E s there are two different 
types of instantons, their actions being different by the 
factor 2. The main contribution to the exponentially 
small subgap DOS is determined by the Gaussian inte- 
gral near the least-action instanton. 

For a planar (quasi-lD) SNS junction with ide- 
ally transparent SN interfaces it is given by (provided 



G N 2/3 < e < 1) 
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where e = (E s - E)/E g , G N = A-kvDL v L z /L x > 1 
is the dimensionless conductance (in units of e 2 /2~kK) 
of the normal part connecting two superconductors, 
E g = 3.12£ T h, E Th = D/Ll is the Thouless energy, 
and 5 = {vV) -1 is the mean level spacing. Here L x 
is the thickness of the N region, assumed to be larger 
than the superconductive coherence length. It is also as- 
sumed that the lateral dimensions L Vl L z are not much 
larger than L x (otherwise, the instanton solution ac- 
quires additional dimcnsion(s) and the exponent 3/2 
changes, cf. 0] for details). The corresponding mean- 
field (MF) expression above the gap reads || 



<p) MF = 3.72 ( 5- 1 v ^. 



(2) 
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Generally, the functional form of Eqs. ([!]), (B is re- 
tained, whereas the coefficients are geometry-dependent 
and can be found from the solution of the standard Us- 
adel equation [|ll] for the specific sample geometry. In 
any case, the total number of states with energies below 
E g is of the order of one. 

In the present Letter we extend our previous re- 
sults in two different directions. Firstly, we derive ex- 
act expression for the DOS in the energy region |e| <C 1, 
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without using anymore the inequality e 3> G ] 
obtained result interpolates smoothly between the semi- 
classical square-root edge (||) and exponential tail ([!]). 
Secondly, we consider the same SNS system allowing 
for non-ideal transparencies at the SN interfaces. The 
result depends upon the relation between dimension- 
less interface conductance Gt and normal conductance 
Gn- As long as Gt > Gj/ 4 , all qualitative features 
of the previous solution are retained, but the value of 
the semiclassical threshold E g and numerical coefficients 
in the expression like (|l]) become dependent upon the 
value of t = Gt/Gn- However, at further decrease 
of interface transparency, Gt <C Gj/ 4 , the DOS be- 
haviour changes dramatically: in the semiclassical re- 
gion E > E g it acquires the inverse-square-root singu- 
larity, (p) MF ~ (— e) -1 / 2 . At smallest |e| this singularity 
smoothens out and crosses over to an exponentially de- 
caying tail of low-energy states. Distinctive feature of 
this tail, as opposed to the situations discussed previ- 
ously, is that the total number of subgap states becomes 
large and grows as G~ 1,2 G)l 8 > 1. We coined this situ- 
ation as "strong tail" , and find exponential asymptotics 
of the DOS in the strong tail region. 

2. Outline of the method. We treat the problem 
within the supersymmetric formalism. The derivation 
of the (T-model functional-integral representation can be 
found in Rcfs. @ [|[ \\J§ . The DOS is given by the in- 
tegral over the supermatrix Q: 



(p(E,r)) = stv(kAQ(r))e- s ^VQ, (3) 

S[Q] = y J dr str [L>(VQ) 2 + AiQ(AE + ir x A)] . 



(4) 



Q is an 8 x 8 matrix operating in Nambu, time-reversal 
and Bose- Fermi (supersymmetry) spaces. Pauli ma- 
trices operating in Nambu and TR spaces are denoted 
Tj and Ui. The matrix k is the third Pauli matrix in 
FB space. A = t z <j z . Integration in (||) runs over the 
manifold Q 2 ~ 1 with the additional constraint 
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Q = CQ 1 C 
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This manifold is parameterized by 8 commuting and 8 
anticommuting variables. It turns out however that only 
4 commuting and 4 anticommuting modes are relevant 
in the vicinity of the quasiclassical gap while contribu- 
tions from all other modes to the DOS cancel. The 
detailed discussion of this fact will be published else- 



where |17| . The reduced parameterization for the com- 
muting part of Q in terms of the 4 variables reads [Q] : 

Qc B = Wz c °sk B + t z smk B (a x cosxb + cr y suixb)] 

x [t z cos 6> b + o- z t x sin B ] , (6) 
Q FF = t z a z cos #p + t x sin 9p . 

The commuting part of the action (with all Grass- 
mann variables being zero) is simplified by introducing 
new variables a = (9 B + k B )/2, /3 = (6 B - k B )/2. Then 
the action (jij) for the normal part (A = 0) takes the 
form 



S[0v,a,0\ = 2S q [9f] - S [a] - S [/3], 



(7) 



S [0] = J dr [D(V6f +4iEcos0] . (8) 

Variation of this action yields the identical Usadel equa- 
tions for Of, a, and [3: 



DV 2 6 + 2iE sin 9 = 0, 



(9) 



with the boundary conditions 9(±L x /2) = n/2. Eq. (Q) 
generally possesses two different solutions 6*i i2 = tt/2 + 
i4>i.2 which coincide (-01, 2(1*) = tpo( r )) just at the thresh- 
old energy E g , and are close to each other in the range 
|e| <C 1 we are interested in. Thus there are 8 possible 
saddle points for the action (^) corresponding to two so- 
lutions of the Usadel equation for each variable 9p, a, (3. 
Rotation over the angle xb € [0,27r) connects some of 
them and produces the whole degenerate family of sad- 
dle points (see Refs. ^, [l7| for details). In the following, 
we will need the function /o(r) which is the normalized 
difference fair) — tpi (r) at E — > E g ; it obeys the linear 
equation 



DV 2 f + 2£ g sinhVo/o = 0. 



(10) 



3. Exact result for the transparent inter- 
face. For energies close to E g we substitute 9 = 
n/2 + iipo + igfo into (|^), expand it in powers of g and 
e and integrate over space using ([K 



S o [0\=S o [*/2 + iil>o] + G 

TTC 2 E g 



G = 



2<S 



-eg- 
2ci 

C2 



(11) 



where we have introduced the constants c„ = 
J(dr/V)fQ n ~ 1 cosh^o- For the quasi-lD geome- 
try, ci = 1.15 and c-i = 0.88. To describe deviation 
of the angles a, (3, 9p from n/2 + iipg we intro- 
duce, analogously to g, three parameters u, v, 



Density of States in a Mesoscopic SNS Junction 



3 



w, respectively. Grassmann variables are intro- 
duced as Q = e -iWe/i e ^W a /2 Ae M a /2 e iW B /2 where 

Q c = e - iW °/ 2 Ae iW °/ 2 is specified in Eqs. (0) and 





jT x CT x (W™) T CT y T : 




as it must satisfy the antiselfconjugate condition W a 
CW^C T = 0. Finally, 



W. 
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Expanding the action in u, v, w and Grassmann vari- 
ables leads to 



S = G 



e(u + v — 2w) 



2w 3 



3 

u + w 



Xn 



v + w 



(12) 



For calculating the DOS we also need an expansion of 
the pre-exponential factor in (^) as well as the Jacobian 
J for the parameterization of the Q-matrix: 



~ / dr str(fcAQ) = + v + 2w), 

8*G 2 
J = 1 u — v I . 



Integrating over Grassmann variables and the 
cyclic angle xb, performing a rescaling (it, v,w) — * 
which excludes G from the inte- 
grand, and changing the variables to I = (u + v)/2, 
m = (u — v) 2 /2, we arrive at the following expression 
for the integral DOS: 



(p) = 



ci(2G) 



-1/3 



Re J dm J dl dw (w + I) 



x (w 2 +2lw+l 2 -m) exp 



l 3 

— — + ew + — - + ml — el 



(13) 



where we introduced the notation e = (2G) 2 / 3 e. 

At this stage we have to choose the contours of in- 
tegration over w and I. The usual convergence require- 
ments for the nonexpanded action (^) enforce the con- 
tour for w (I) go along the imaginary (real) axis at large 
values of w (I). However since the main contribution 
to the DOS comes from the expression (Ot) determined 




Fig. 1. Possible contours for integration over w and 
The proper choice is Ci for w and C3 for I. 



by small w and I, these contours should be properly de- 
formed to achieve convergence of ( |l3| ) . The integral ( [l3"|) 
converges if the contour for I runs to infinity in the dark 
regions in Fig. 1 and otherwise for w. Therefore, we 
should choose the contour C\ for w (see Fig. 1), whereas 
for I there are two possibilities: Ci and C3. The correct 
choice is dictated by positivity of the DOS, which imply 
the contour C3 for /. Integration in ( |l3| ) is straightfor- 
ward although rather cumbersome and leads to the final 
expression for the DOS: 



(P) 



2-7TC1 

(2G) 1 /3 < 5 



- £ Ai 2 ( £ ) + [Ai'( e )] 2 
Ai(e) 



+ 



dy Ai(y) 



(14) 



where Ai(e) is the Airy function. Asymptotic behaviour 
of the calculated DOS at e 3> 1 coincides with the result 
(|l]) of the single-instanton approximation j|, see Fig. 2. 

The functional dependence (]lj) coincides with the 
RMT prediction for the spectrum edge in the orthogonal 
ensemble p2| . It is not surprising because the random 
matrix theory is known to be equivalent to the 0D a- 
model In our case the problem became effectively 
0D when we fixed the coordinate dependence /o(r) for 
the parameters of Q near E g . 

Breaking the time-reversal symmetry drives the sys- 
tem to the unitary universality class. The correspond- 
ing RMT result @ can be obtained from Eq. @ by 
dropping the last integral term. This result can be eas- 
ily derived by the cr-model analysis in the following way. 
Strong magnetic field imposes an additional constraint 
on the Q-matrix. As a result the mode associated with 
the variable m acquires a mass, so, instead of integrat- 
ing over it we set to = 0. One of the Grassman modes is 
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Fig. 2. Exact dependence ( p^ ) of the DOS on the di- 
mensionless energy e (solid line). Single-instanton ap- 
proximation (^) (dotted line). Semiclassical result (g) 
for the DOS (dashed line). 



also frozen out giving the pre-exponent (w + I) 2 in the 
integral ( [13|) . Finally, the expression for DoS coincides 
with Eq. (|14| ) but without the last term. 

4. Finite transparency of the NS interface. 

Now we turn to the analysis of the sub-gap structure of 
a quasi-lD SNS contact with finite conductance Gt of 
the NS boundary. The role of the interface is described 
by the dimensionless parameter t = Gt/Gn- For t 3> 1, 
the interface is transparent and the result (|l4|) apply. 
In what follows we will consider the case t <C 1. The ef- 
fect of finite transparency is described by the additional 
boundary term [fL4j in the action 



1G 



str (Q L Q R ) 



(15) 



where Q L,R are the Q-matrices at both sides of the in- 
terface. Eq. ( |l5| ) is the the first term in the expansion 
of the general boundary action JlJ, |J [ll| in the small 
transparency T <C 1 of conductive channel and leads 
to the Kupriyanov - Lukichev |l^] boundary conditions. 
In the diffusive regime (I <C L x ) at t <C 1, we have 
r ~ tl/L x that justifies the use of Eq. (|TJ). The com- 
muting part of the action can still be written in the 
form (^) with the additional term in 5*0 : 



So[0] = 



J dx [D{9'f +UE cos ( 



-W 2 



Gi 



■ sin 



~2~ 



(16) 



In the limit t -C 1, the Usadel equation has almost 
space homogeneous solutions, which allows to use an ex- 
pansion tp = A+B[l — A(x/L x ) 2 ] for them. Substituting 



this ansatz into the action (jig), and minimizing over B 
we obtain the action in terms of P = e A : 



So(P) = % 



2t t 2 9 

(s-t)P + — P 2 

v ' P 24 



(17) 



where s = E / '-Exh- Here we keep only leading terms and 
substitute all s except the first one by t. After variation 
we find the cubic saddle point equation for P: 



1 



2 

P 2 



f 

12' 



-P. 



(18) 



The maximum of the RHS achieved at P = (48/i) 1 / 3 
determines the position of the mean-field gap: s g = 
t + 0{t b / 3 ) and hence E g — GtS/^tt. Depending on 
the deviation from the threshold, e = (E g — E)/E g — 
(s g — s)/s g , there are two regimes for Eq. ([l8|). 

Weak tail. If |e| < t 2 ^ the two solutions of @ 
are close to each other and can be seek in the form 
P = Pq + SP. Expanding the action in powers of 5P we 
get 



S (P) = So(Po) + ^ 



z5P + —,{5Pf 

■Ml 



(19) 



This equation closely resembles its counterpart (pd| ) for 
the transparent interface. As mentioned in Ref. Q7this 
form of the expansion of the action over small devia- 
tions near the threshold solution inevitably leads to the 
instanton action scaling as e 3 / 2 . In fact, there is full 
equivalence |i~7j between the DOS for the transparent 
NS interface given by Eq. ([IJ) and the DOS in the limit 
lei < i 2/3 < 1. The latter can be obtained from the 
former by redefinition of the constants c\^- For a ID 
planar contact they appear to be ci = Po/2, ci — 6/Po, 
E g / i?Th = t. 

In particular, above the threshold, at e < 0, one 
encounters the mean-field square-root singularity 



(P) 



MF 



4 Ql/6 



(20) 



The instanton action becomes S = Sq(Pi) — Sq(P2) = 
(2/3)6 1 / 6 GN^ 1 ^ 3 e 3 ^ 2 , and the single-instanton asymp- 
totics of the DOS tail reads 



1 / 7T6 1 / 6 

{P) ~ 6 V 2G N i 5 / 3 v^ 



exp 



■~6 1 '*G N t L W ) . (21) 



Strong tail. In the opposite limit, t 2 / 3 -C |e| <C 1, 
the difference between the two solutions for Eq. ( p"8| ) 
is large but expansion (|l^) is still valid (gradients of 
■01.2 are small provided e <C 1). The roots P% t 2 can 
be found neglecting either the second or the third term 
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in Eq. @: P 1 = yffle, P 2 = 12e/t, with P 2 > J\. 
Above the threshold that gives the mverse-square-root 
singularity in the semiclassical DOS: 



(p)mf = - Re 



dr cos6> = ilm(P- P -1 ) = \ 






Below E g one obtains for the instanton action S = 
—So(P2) = (3/4)Gn£ 2 which determines the one-in- 
stanton asymptotics of the subgap DOS. The pre-ex- 
ponent can be calculated by generalizing the method of 
Ref. |Q. Introducing the deviation parameter q accord- 
ing to a = 7r/2 + ilogP2 + iq/V2 and expanding the 
action in powers of q and the corresponding Grassmann 
pair we obtain for the action and the pre-exponential 
factor in Eq. (0): 



S — -Gn£ 2 
o 



2-q 2 



2\/2. 



~ j dr str(fcAQ) = 



i-4= 



V2 



The measure of 
2yJWft{2ey/ i dqdC,d£ > . 
we finally obtain 



integration is T>Q 
Inserting these into Eq. 



yexp 



(23) 



5. Discussion. We have considered the integral 
density of states in a coherent diffusive SNS junction 
with arbitrary transparency of the SN interface. For 
the ideal interface (Gt 3> Gn) we managed to go be- 
yond the single-instanton analysis Q and derived the 
exact result ( fl4| ) valid as long as \E — E s \ -C P g . This 
expression uniquely describes the semiclassical square- 
root DOS (||) above the Thouless gap E s , the far subgap 
tail (|l|), and the crossover region e ~ G^ 2 ^ 3 between the 
two asymptotics. The functional form of this result co- 
incides with the prediction of the RMT. 

As the SN interface become less transparent, Gt <§; 
Gn, the situation changes. At Gt Gj!/ 4 these changes 
are only quantitative: the position of the quasiclassical 
gap is shifted to E g = (G T /G N )E Th , but the DOS both 
above [Eq. (po|)1 and below [Eq. (fjlj)] the gap has the 
same dependence on the deviation e from E g , with the 
coefficients becoming dependent on Gt- In this limit, 
the very far part of the tail [at e » (Gt/Gn) 2 / 3 ] ex- 
hibits another e-dependence (f23"[), but the correspond- 
ing DOS is exponentially small. Therefore, in the limit 

1 /4 

Gt ^> G N the total number of the subgap states is of 
the order of 1 and is independent on Gt- We refer to 
this case as weak tail. 





e < 


£ > 


£ fluct 


t » 1 


(i 


) 




r -2/3 


\s\ < t 2 ^ < 1 


© 




G N 2/3 <- 2 / 9 


i 2/3 < |e| < 1 


(1 









Table I. References to the asymptotic formulas for the 
DOS above (e < 0) and below (s > 0) the gap, and the 
width |e|fl uct of the fluctuation region for the regimes of 
the transparent interface (t S> 1), weak (|e| <C t 2 ^ <C 1) 
and strong (t 2/3 < |e| < 1) tails. 



As the interface becomes less transparent, the region 
of applicability of the weak tail shrinks and finally dis- 
appears at Gt ~ Cj/ 4 . For even lower G T < G?/ 4 , 
the difference between the case of the transparent inter- 
face becomes qualitative: the DOS above E g acquires 
an inverse square-root dependence (p^), while the sub- 
gap DOS follows d25|). In this regime the total number 
of the subgap states is proportional to G^^Gj/ 8 3> 1 
and grows with decreasing Gt in contrast to all previ- 
ous cases when this number is of the order of 1. This 
indicates that at Gt ~ Gj/ 4 the universality class of the 
problem is changed. At Gt <C Gj/ 4 it is no longer equiv- 
alent to the spectral edge of the Wigner-Dyson random 
matrix ensembles. 

The asymptotic results for the DOS above and below 
the gap, as well as the width of the fluctuation region 
near E g are summarized in Table I for the three regions 
considered. 
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